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Page 30, caption, for ease of use please note the following conversions: 

0.05 N sec/m^ = 50 cP 
0.60 cm^/sec = 60 cSt 

Page 30, table III, Xa formula should read as follows: 
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Page 30, table III, Xp formula should read as follows: 
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Page 32, table III, Vs formula should read as follows: 

0. 1 047 x([ + mg)xnpX(X -Xp)/mg 


Note: We gratefully acknowledge C.E. Douglas for bringing these corrections to our attention. 



SUMMARY 


A simple method for predicting the power loss of a steel spur gearset of arbitrary 
geometry supported by ball bearings was developed. The method algebraically accounts 
for losses due to gear sliding, rolling traction, and windage and incorporates an expres- 
sion for ball-bearing power loss to obtain an accurate estimate of spur-gear-system 
efficiency at part-load as well as full-load conditions. The analysis was compared with 
data from an experimental investigation. Theoretical predictions generally agreed with 
the test data except at low oil flow rates, where the analysis underestimated the gear- 
system efficiency. The theoretical sensitivity of power loss to changes in speed and 
load was close to that of the experimental data. A theoretical comparison of the contri- 
bution of the individual component losses relative to the total system losses indicated 
that losses due to rolling traction, windage, and support bearings are significant and 
should be accounted for in the calculation of power loss. 


INTRODUCTION 

Investigations of gear power loss are numerous, as evidenced by the research re- 
ported in references 1 to 5. From a practical standpoint, few, if any, of these analyti- 
cal methods were intended to estimate gear efficiency accurately at less than full-load 
conditions. It is generally true that a gearset that is fully loaded will provide higher 
efficiency than one that is partially loaded. Since most machines are sized for occasion- 
al overloads, they operate at power levels substantially lower than maximum for most 
of their lives. With today's emphasis on reducing energy consumption of machinery, it 
is important to be able to accurately determine power losses at less than full-load con- 
ditions. 

The basic limitation of most of the existing methods of predicting gear efficiency is 
that they rely solely on the selection of a coefficient of friction to determine power loss. 
There are various methods available for calculating this friction coefficient, each result- 
ing in a different value. When part -load gear efficiency is calculated, the uncertainty in 
the selection of this coefficient may overshadow the variation in power loss at the part- 
load condition. 

In addition, the friction coefficient, which can account for tooth sliding losses, can- 
not be used to calculate the rolling (or pumping) loss also present in a gear mesh. The 



rolling loss is caused by hydrodynamic forces on the gear teeth. This loss, as well as 
gear -windage and support -bearing losses, can become a significant portion of the total 
mesh loss, particularly at light loads and high speeds. 

The spur-gear efficiency analysis of Y. P. Chiu (ref. 5) includes the effects of both 
sliding and rolling in the gear mesh. Expressions for rolling and sliding losses as ap- 
plied to spur -gear geometry are based on roller test data. Instantaneous values of slid- 
ing and rolling power loss are integrated over the path of contact of the gear and then 
averaged to obtain an average power loss. The method of reference 5 is perhaps the 
most thorough to date. However, it does not include gear -windage loss, support- 
bearing losses, or the effects of tooth load sharing across the mesh cycle. The support- 
bearing and rolling losses, as shown later in this report, largely dictate the part-load 
gear-system efficiency. Also, the expression for the coefficient of friction used in this 
work predicts an abnormally high loss when compared with test data. 

The objectives of the analysis reported herein are (1) to provide a general method 
for calculating the full- and part-load efficiency of a spur -gear system of arbitrary 
geometry which considers losses due to gear sliding, rolling, and windage and rolling- 
element support -bearing losses and (2) to compare the predictions from this analysis 
and the analysis of reference 5 with spur-gear performance test data of reference 6. 
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semimajor width of contact ellipse, m (in. ) 
semiminor width of contact ellipse, m (in. ) 
dimensionless moment coefficient 

dimensionless moment coefficient for disk of zero thickness 

support -bearing basic static capacity, N (Ibf) 

constants of proportionality; see table II 

pitch circle diameter, m (in. ) 

tip diameter, m(in.) 

bearing pitch diameter, m (in. ) 

base circle diameter, m (in. ) 

modulus of elasticity, N/m^ (Ibf/in^) 


equivalent modulus of elasticity, 
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normal applied load, N (Ibf) 

rolling traction force, N (Ibf) 

sliding force, N (Ibf) 

static equivalent bearing load, N (Ibf) 

combined radial and tangential bearing load, N (Ibf) 

face width of tooth, m (in. ) 

ball-bearing lubrication factor 

coefficient of friction 

dimensionless material parameter, £' 0 - 

dimensionless film thickness (eq. (4)) 

isothermal central film thickness, m (in. ) 

thermally corrected film thickness used in rolling-traction equation 
(eq. (5)), m (in.) 

constant defined in appendix C 

lubricant coefficient of thermal conductivity, W/m K (Btu/hr ft °F) 
ellipticity parameter, a/b 

path of contact distances defined in appendix A, m (in. ) 
bearing friction torque, N-m (in-lbf) 

load -dependent part of bearing friction torque, N-m (in-lbf) 

viscous part of bearing friction torque, N-m (in-lbf) 

gear ratio, Ng/Np 

number of gear teeth 

efficiency, percent 

rotational speed, rpm 

power loss 

total power loss due to rolling-element support bearings, kW (hp) 

power loss due to rolling traction, kW (hp) 

power loss due to tooth sliding, KW (hp) 

power loss due to windage, kW (hp) 

diametral pitch 
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2 2 

lubricant pressure, N/m (Ibf/in ) 
base pitch, m (in. ) 

dimensionless thermal loading factor, 6/K^ 

pitch circle radius or radius in general, m (in. ) 

equivalent rolling radius, m (in. ) 

radius at point of contact, m (in.) 

effective radius in direction of rolling, m (in. ) 

rotating disk Reynolds number, C^^pcuR /u 

axial clearance between gear and housing wall, m (in. ) 

pinion torque, N-m (in-lbf) 

disk thickness, m (in. ) 

dimensionless speed parameter, ut?q/E'R^ 

average rolling velocity, (V + V )/2, m/sec (in/sec) 

o 

surface velocity, m/sec (in/sec) 
sliding velocity, V - V , m/sec (in/sec) 

O 

rolling velocity, V + V , m/sec (in/sec) 

O 

dimensionless load parameter, Fjj/E’R^ 
gear contact normal load, N (Ibf) 
maximum normal gear tooth load, N (Ibf) 
transmitted load, N (Ibf) 

path of contact distances defined in appendix A, m (in.) 

exponent in Mp equation based on bearing geometry (eq. (15)) 

factor in equation based on bearing geometry (eq. (15)) 

2 2 

pressure-viscosity coefficient of lubricant, m /N (in /Ibf) 
temperature, K (°F) 

Poisson's ratio 

temperature -viscosity coefficient, 1/K (1/°F) 

Reynolds number exponent 


4 



e 

X 




M 


air 


eq 


P 

^air 

^eq 

U) 

Subscripts: 


gear tooth pressure angle, deg 

dimensionless ratio of film thickness to composite surface roughness 

-3 2 2 

lubricant absolute viscosity, 10" N sec/m (cP) (Ibf sec/in ) 

-3 2 2 

absolute viscosity of air within gearbox, 10“ N sec/m (cP) (Ibf sec/in ) 

-3 2 2 

equivalent air-oil absolute viscosity, 10" N sec/m (cP) (Ibf sec/in ) 

-2 2 2 

lubricant kinematic viscosity, 10" cm /sec (cS) (ft /sec) 

3 3 

lubricant density, kg/m (Ibm/ft ) 

3 3 

density of air within gearbox, kg/m (Ibm/ft ) 

3 3 

equivalent density of air-oil mixture, kg/m (Ibm/ft ) 
thermal reduction factor 
angular velocity, rad/sec 
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bearing 

gear 

input 

pinion 

rolling 

sliding 

total 

ambient 


conditions 


Superscripts: 

(■“) average value 

(") simplified value 


ANALYSIS 

The most significant sources of power loss in a spur -gear system include the gear 
mesh, gear windage, and support bearings. In this analysis, is assumed that the gears 
are made of steel and that they are jet lubricated. Furthermore, it is assumed that the 
gear does not come into contact with oil in the sump, so that there are no oil churning 
losses. Secondary losses such as momentum transfer due to oil impinging on the gear 
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teeth or noise generation are not treated. The windage and bearing losses can be cal- 
culated in a straightforward manner with approximate expressions. The mesh losses 
are more complex and are analyzed in detail. 


Gear Mesh Losses 

The mesh losses consist of a sliding frictional component and a hydrodynamic roll- 
ing component. The hydrodynamic rolling (or pumping) loss is the power required to 
entrain and compress the lubricant to form a pressurized oil film which separates the 
gear teeth (ref. 7). It is shown later that at light loads the rolling traction loss is a 
major portion of the system loss. 

The calculation of sliding and rolling power loss is based on data from roller test 
machines, where friction due to sliding and rolling have been measured under various 
operating conditions. The gear contact progressing through the mesh cycle can be 
modeled as a constantly changing roller contact whose size, speed, and load can be cal- 
culated from involute gear geometry at the given operating conditions. The gear mesh 
cycle is then approximated by a large number of discrete elastohydrodynamic (EHD) con- 
tacts. The instantaneous rolling and sliding forces can be calculated for each contact, 
and the results integrated over the path of contact for the complete mesh cycle to obtain 
an average gear contact power loss. 

The path of contact, the abscissa in figure 1, provides a convenient coordinate sys- 
tem for calculating the gear mesh losses. Figure 1 shows the assumed instantaneous 
gear tooth loading as three sequential pairs of gear teeth come into contact. The coordi- 
nates called out on the path of contact represent points where teeth enter or leave a 
mesh cycle. The path of contact is actually a line in space of fixed inclination (the pres- 
sure angle) mapped out by the instantaneous contact point between two mating gear teeth. 
During the mesh cycle, the load transmitted between the gears is normally carried by 
either one or two teeth at any one time. The average number of teeth in engagement is 
generally referred to as the contact ratio. Figure 1 illustrates the tooth loading pattern 
used in this analysis. Two teeth share the load for most of the path of contact. This 
gear mesh loading is based strictly on involute gear geometry with no profile modifica- 
tion. Dynamic loading is not included here, but its effect on efficiency is judged to be 
small. 

Figure 2 shows the variations in the gear mesh velocities and the equivalent rolling 
radius as the teeth move through the mesh cycle for the spur-gear geometry used in ref- 
erence 6 (see table I). The data of reference 6 are used later for comparison with the 
present analysis. The variations in these parameters show that the conditions in a spur- 
gear mesh are constantly changing, and, thus, the instantaneous losses are changing as 
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well. Appendix A gives the equations used to calculate points along the path of contact, 
gear loads, sliding velocities, rolling velocities, and values of the equivalent rolling 
radius. 

Sliding force . - The instantaneous frictional force due to sliding of two gear teeth 
against each other is 


Fg(X) =/(X)w(X) (1) 

The selection of the coefficient of friction can significantly affect the system losses, 
since gear tooth sliding is a major source of power loss in a loaded gearset, as shown 
later. For this analysis, the coefficient of friction developed by Benedict and Kelley 
(ref. 8) is used: 


CjW(X) 

/(X) = 0.0127 log 'L 

/io[Vs(X)][VT(X)]2 


( 2 ) 


where 


Cj^ = 29. 66 (SI units) 

= 45.94 (U. S. customary units) 

(A summary of the constants used in all equations is presented in table n. ) This coeffi- 
cient is probably the most accurate expression available at present for use when the true 
gear surface temperature is not known (ref. 7). There is no universal agreement as to 
the proper form of the friction-coefficient equation, particularly with regard to the ef- 
fect of inlet oil viscosity. At low values of A (ratio of minimum EHD film thickness to 
composite surface roughness), a significant portion of the friction is due to asperity con- 
tact. An increase in viscosity decreases the number of asperity contacts, and, thus, 
the friction coefficient decreases, as indicated by the Benedict and Kelley expression. 

If, on the other hand, full-film lubrication exists, an increase in viscosity increases the 
friction coefficient. As a practical matter, very few gear applications would be ex- 
pected to have full-film lubrication. 

Benedict and Kelley conducted their tests on rollers with an initial rms surface 
roughness of 0. 4 micrometer (16 pin.). The gears of reference 6 (cited in this report 
for comparison with theory) were precision hobbed, a process that results in an rms 
surface roughness of approximately 0. 4 micrometer (16 pin.). Since the operating 
speeds and values of A in reference 6 are similar to those in the tests done by Benedict 
and Kelley, the friction equation (eq. (2)) should be applicable. Strictly speaking, equa- 
tion (2) is valid when A is less than approximately 2, which is the range of most gear 
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applications. Extending its use to systems with higher values of X should not introduce 
serious errors. 

The variation in coefficient of friction as determined by equation (2) for the gear 
geometry of reference 6 is shown in figure 2. The coefficient of friction varies along 
the path of contact because of the variation of the sliding velocity. Equation (2) does not 
predict accurate values of the friction coefficient at low sliding velocities (near the pitch 
point). An upper limit of f =Q.\ was set in the computer program. This limit should 
not introduce serious errors in the power loss calculation since the sliding power loss 
is relatively low at this point and of short duration. 

Rolling force . - The calculation of rolling traction is common to bearing analysis 
(ref. 9), but is generally overlooked in calculations of gear power loss. The resistance 
a ball encounters when rolling over a lubricated bearing race is due to the hydrodynamic 
pressure forces that build the EHD film in the ball-raceway contact. Similar phenomena 
occur as two rollers roll over each other or as two gear teeth move across one another. 
The buildup of pressure on the rolling surfaces resists the motion of the two bodies and, 
therefore, absorbs power. 

Experimental work conducted by Crook (ref. 10) on rolling traction in roller con- 
tacts can be used to calculate the rolling losses in a gear mesh. The rolling traction 
force can be expressed as 




(3) 


where 

C 2 = 9 . 0 x 10 "^ (SI units) 

= 1.3x10^ (U. S. customary units) 

In this investigation, it is assumed that the gear contact operates principally in the 
viscous -elastic lubrication regime. The gear contact film thickness is calculated by the 
method of Hamrock and Dowson (ref. 11); 

H„ =2.69 l/’- (1 - 0. 61 e"®- "^^^) (4) 

Rx 

(See appendix B for application of this equation to the present analysis. ) At high pitch 
line velocities, this equation generally overestimates the actual film thickness. The 
reason for this is that thermal and starvation effects, which are not included in the equa- 
tion, limit the buildup of the lubricating film. A thermal reduction factor developed 
in reference 12 (see fig. 3) is used to limit h at high speeds as follows: 

hj^ = h^^ (5) 
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Therefore the instantaneous rolling traction force for a gear rriesh becomes 

Fj^(X) = C^hiXW^iX)-/ (6) 

Sliding and rolling power loss . - The instantaneous sliding and rolling power loss 
can be expressed as 


Pg(X) = C3Vg(X)Fs(X) ' 
Pj^(X) = C3 V^(X)Fj^(X) 


( 7 ) 


where 

C3 = 10"^ (SI units) 

= 1. 515x10^ (U.S. customary units) 


The average power loss over one mesh cycle is obtained by integrating Pg(X) and 
Pj^(X) over the path of contact and then dividing by the length of the path of contact. 

Since the functions Pg(X) and Pj^(X) do not lend themselves to a closed-form analytical 
solution, a Simpson's rule numerical integration technique is applied. Figure 1, mesh 2, 
illustrates the implementation of the power loss integration. To obtain the total power 
lost during the mesh 2 cycle, we must include not only the losses incurred in mesh 2 but 
also the losses in mesh 1 from point X3 to X^ and the losses in mesh 3 from point Xj^ 
to Xg, since the load is being shared by the teeth just entering and leaving engagement 
during these periods of time. Therefore 


Pc + P 


R 


X. -X, 



Pj^(X)]dX 


+ 



[Pg(X)+ Pj^(X)]dX + 2 



[Pg(X) + Pj^(X)]dX 


( 8 ) 


In the numerical integration of these expressions, the sections of the path of contact 
between Xj^ and X2 and between X3 and X^ are each divided into 200 equal steps. 
The path of contact between Xg and X3 is divided into 40 equal steps. Calculation of 
the average power loss at various step sizes produces convergence of the solution even 
with very few steps. Thus the use of the step size indicated yields an accurate integra- 
tion of the functions. 
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Windage Losses 


Gear windage is generally a smaller portion of the overall gear-system loss at low 
to moderate speeds than at high speeds. (Windage is defined here to be the power re- 
quired to rotate the pinion and gear in the air-oil atmosphere present within the gearbox 
It does not include churning losses, since the gears are assumed not to come in contact 
with an oil reservoir, nor does it include oil fling-off losses. ) Windage increases rapid 
ly with speed. From work done on turbine rotor wheels of references 13 and 14 (see 
appendix C), the following expression for windage power loss can be written: 


W 


l + 2.3lV'-V- W-2 

R/ 


(9) 


where 

= 2.04x10'® (SI units) 

= 2. 70x10”^® (U. S. customary units) 


The environment within a gearbox is that of an air-oil mixture. Thus the density and 
viscosity in equation (9) must represent this mixture. In reference 15, an equivalent 
air -oil mixture density is calculated to simulate a typical environment within a helicop- 
ter gearbox, an environment similar to that in most jet-lubricated gearboxes. The 
equivalent density is calculated for a mixture having a nominal composition of 34.25 
parts air to 1 part oil (an oil/ air ratio of about 3 percent). For this analysis, both the 
density and the viscosity are calculated to represent this mixture. Thus 


eq 


eq 


(l)p + (34.25 )p^.^ 
35.25 

(Dm + (34.25)p^.^ 
35.25 


( 10 ) 


For this analysis, all properties in equation (10) except oil viscosity are assumed con- 
stant. Air density and viscosity at 339 K (150° F) are assumed to be 

p_,„ = 1.042 kg/m® (0.065 Ibm/ft®) 


^ = 2x10'® N sec/m® (2.94x10"® Ibf sec/in®) 

cllx 

Assuming a specific gravity of 0. 9 for the oil yields 


10 



p = 899. 3 kg/m^ (56. 2 Ibm/ft^) 


Use of these values results in 


pgq = 26. 52 kg/m^ (1. 66 Ibm/ft^) ' 


p = 0. 028 p + Cr 
eq ' 5 


where 


( 11 ) 


= 0. 019 (SI units) 

= 2. 86x10”^ (U. S. customary units) 


Combining equation (11) with equation (9) and using gear notation give the windage power 
loss for the pinion and gear: 


^\\T „-C^(l + 2.3 — 


W,g -'-6 




, 2 . 8 


VV"gy 


®(0.028 p + 


Pw, p = Cg (1 + 2. 3 ^ j n2- 8 r4. ^28 p + C^)^- ^ 


( 12 ) 




where 


Cg = 2. 82x10"^ (SI units) 

= 4.05x10"^^ (U. S. customary units) 


Rolling-Element Support-Bearing Losses 

Bearing power loss is a much investigated and diverse subject in itself. However, 
in this investigation, we are interested in obtaining an engineering estimate of the bear- 
ing power loss without the need to resort to a full computer analysis. Accordingly the 
simple method described by Harris (ref. 9) is used. Of course, other methods or data 
could be used if available. It is assumed here that the pinion and gear are straddle 
mounted and centered on the shaft with identical-size deep-groove ball bearings. 
Furthermore the pitch diameter is assumed to be the arithmetic average of the bearing 
bore diameter and the outside diameter. These dimensions can be readily found in a 
bearing catalog along with the bearing static capacity. From reference 9, the friction 
torque for one bearing is made up of a load-dependent term and a viscous term: 
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M = Mj^ + My 


(13) 


The load-dependent term is 


Ml =zi 


^ST 




(14) 


For a deep-groove bearing, z = 0. 0009 and y = 0. 55. In the case of equally spaced, 
straddle -mounted spur gears, no significant axial forces are developed, so that 
= Fg^, and equation (14) becomes 


Ml = 0. 0009 


.1. 55 
ST 

,0.55 


D 


m 


(15) 


The viscous term My is given as, 

for m > 2000 (SI units) (0.0215 (U.S. customary units)), 


M 


V = 1.42xl0"^f C„ 

V o' ' m 7 


(16a) 


for t^n s 2000 (SI units) (0.0215 (U.S. customary units)), 

M^ = 3.492xl0”^f D^Cn 

V 0 111 o 


(16b) 


where 


= 6894 (SI units) 

= 2051 (U.S. customary units) 

Cg = 6894 (SI units) 

= 1.0 (U.S. customary units) 

and f^ is a term used to represent the amount of lubricant in the contact. An f^ value 
of 2 is selected; it corresponds to that for a partially flooded contact, such as an oil 
bath or jet-lubricated bearing, an arrangement halfway between a mist and flooded oil 
lubrication. The torque loss is calculated for both the pinion and the gear bearings and 
doubled, since there are two bearings on each shaft. The torque loss is then converted 
to a power loss by multiplying by the appropriate bearing speed: 
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where 


Cg = 1.05x10'^ (SI units) 

= 1. 59x10”^ (U. S. customary units) 


RESULTS AND DISCUSSION 


Simplified Equations for Sliding and Rolling Power Loss 


Figure 2 shows the calculated variations of sliding and rolling power loss along the 
path of contact for the gearset of reference 6 at the full-power test point. It must be 
noted that the losses include contributions from two sets of teeth when the transmitted 
load is shared. The simple shapes of these functions suggest that algebraic expressions 
can be developed to accomplish the integration, so that the equations for sliding and 
rolling loss are simplified. These simplified expressions are 




lPs(«i) + Ps(<2)]«3 + 


IPs(«4)j«5) 


(18) 


p [PR<V-^PR<yl^3-^tPRgp)l«5 ( 1 , 

The results from the simplified equations are within 1 percent of the numerically inte- 
grated results for a large number of cases. 


Comparison of Analysis with Data 

The accuracy of the present analysis was validated by comparisons with experi- 
mental spur-gear performance data. Part -load gear efficiency data are very limited 
largely because of the difficultly in accurately measuring small power losses. In the 
study of reference 6, power loss was measured by the square-loop test method with 
speed, load, oil viscosity, gear width, location of oil jet, and lubricant flow rate as test 
variables. Although the present analysis does not account for oil flow rate or jet loca- 
tion, the effects of speed, load, and face width are included and were compared with 
test data. The support -bearing loss was also measured in the study of reference 6, so 



that the gear power loss could be separated from the bearing loss. The gear geometry 
and operating conditions used in this comparison are listed in table I. 

The appropriate oil viscosity to be used in the calculations requires an estimate of 
the bulk gear tooth temperature at the surface, since calculations of the sliding friction 
coefficient and the EHD film thickness are based on this temperature. Unfortunately the 
bulk gear tooth surface temperature is rarely known, so that an estimate must be made. 
A reasonable approximation can be made by assuming that the bulk tooth surface tem- 
perature (not the same as the instantaneous tooth temperature along the path of contact) 
is equal to the temperature of the oil exiting from the mesh. In the tests of reference 6, 
only the oil supply temperature of 333 K (140° F) was given. This temperature was used 
to make an initial estimate of the mesh power loss. The initial calculations indicated 
that the maximum gear losses are approximately 0. 37 kilowatt (0. 5 hp). If the entire 
0. 37 kilowatt (0. 5 hp) were transmitted to the oil at the lowest flow rate tested, 1. 9 
liters per minute (0. 5 gal/min), the maximum mesh exit oil temperature would be 340 K 
(153° F). Therefore, a nominal bulk gear surface temperature of 339 K (150° F) was 

selected, and the approximate absolute oil viscosity to be used in the gear power loss 
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calculations was then 0. 05 N sec/m . Since the temperature rise across the bearing is 
small, the calculations of bearing power loss use the absolute viscosity at the supply oil 
jet temperature of 333 K (140° F), 0.054 N sec/m^. An example calculation for the 
maximum -power test point of reference 6 is given in table III. 

In figure 4 the calculated gear losses (excluding bearing losses) are compared with 
the test data of reference 6 for a gear width of 0. 04 meter (1. 563 in. ) and for various 
lubricant flow rates, oil jet positions, and pinion speeds as a function of pinion torque. 
Also shown is the gear power loss as calculated by the method of reference 5. Oil flow 
rate was held at either 1. 9 or 11.4 liters per minute (0. 5 or 3 gal/min), and the oil sup- 
ply jet was positioned so that the oil impinged on either the inlet or the outlet side of the 
mesh. In all cases, the measured power loss was greatest at the higher flow rate with 
the oil jet positioned at the inlet to the mesh. The current analysis predicts slightly 
higher power losses than measured, but predicts with reasonable accuracy the sensitiv- 
ity of gear mesh power loss to speed and load. The analysis of reference 5 predicted 
significantly greater power loss than the current analysis. The reason for the differ- 
ence between the current analysis and the prediction from reference 5 lies, principally, 
in the friction-coefficient expressions used. In figure 4 as speed is increased, the 
power loss becomes more sensitive to changes in the applied loads. This change in sen- 
sitivity can be understood if the power loss is considered as a loss of torque at the drive 
pinion shaft. At each operating speed, there is some variation in torque loss as a func- 
tion of load, but the variation is small. The power loss is the product of torque loss and 
pinion speed. The pinion speed amplifies the torque loss when expressed as power. 

Most of the change in sensitivity of power loss with speed shown in figure 4 is due to 
amplification of the torque loss by the operating speed. 
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In figure 5, the measured bearing losses of reference 6 are compared with pre- 
dicted results. The calculated bearing losses are approximately 25 percent low at each 
speed, but the load sensitivity is accurate. 

In figure 6, the gear and bearing losses are combined and shown as overall gearbox 
efficiency as a function of pinion torque. The calculated efficiency is within 0. 3 per- 
centage point of the measured results at an oil flow rate of 11. 4 liters per minute 
(3 gal/min) with oil supplied into the mesh. When the analytical results are compared 
with the experimental data for the other lubricant flow rates and oil jet positions, the 
largest difference occurs at low torque and high speed with the most efficient lubrication 
method, namely, a flow rate of 1. 9 liters per minute (0. 5 gal/min) and the oil supplied 
to the exit of the mesh. At this condition, the analysis predicts an efficiency that is 0. 6 
percentage point lower than the measured value. The predicted trends with speed and 
torque generally show good agreement with the data. 

In the tests of reference 6, data were also taken on a gear of larger face width but 
of the same basic geometry (see table I). Figure 7 shows the predicted system effi- 
ciency for the wider gear. The part-load efficiency for this gear drops more substan- 
tially at higher speeds than the narrower gear. The analysis correctly predicts this 
lower efficiency. As in the previous case, the calculated values agree closely with the 
test data for the least efficient method of lubrication. The test data at 250 rpm for the 
wider gear show more variation between the types of lubrication. For this geometry, 
the greatest deviation between the analytical and measured data is 0. 7 percentage point, 
at 2000 rpm and 67. 8-newton-meter (50-ft-lbf) pinion torque. When the analysis is 
compared with the data for 11. 4 liters per minute (3 gal/min), the maximum difference 
is 0. 5 percentage point at 250 rpm and 67. 8 newton-meters (50 ft-lbf). Once again the 
analysis shows remarkably good agreement in speed and load variation with the test 
data. 

In the event that full-load gear -system data become available either empirically or 
through another analytical technique, the present analysis could be calibrated to match 
these data at the full-load efficiency point. The adjusted analysis could then be used 
with a high degree of confidence to predict off-speed and off-load design-point perfor- 
mance. 


Breakdown of Gear -System Losses 

With the gear loss prediction method presented in this study, the theoretical con- 
tributions of the individual component losses to the total gear mesh loss can be readily 
determined at any speed and applied load. An example of this is given in figure 8, which 
shows a theoretical percentage breakdown of the various components of gear-system 
power loss for the 0. 04-meter-(l. 563-in. -) wide test gear as a function of applied load 
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for three speeds. At low speeds, the sliding loss accounts for most of the system loss. 
At higher speeds, the sliding loss becomes less important, and the rolling loss and the 
pinion bearing losses become significant. Gear and pinion windage are often neglected. 
At low speeds, this is justified. However, at higher speeds such as the value 2000 rpm 
used in this example, windage contributes as much as 10 percent of the total power loss 
and should not be neglected. At low torque, the sliding loss is low, since this loss is a 
direct function of load. The rolling loss is proportional to film thickness, which is 
relatively insensitive to load. Therefore rolling loss dominates at low torque. 

Figure 8(c) shows the error in using the sliding loss alone, that is, using just a 
gear coefficient of friction, to predict gear efficiency as suggested by several investi- 
gators (refs. 1 to 3). At the full load of 271 newton-meters (200 ft-lbf), the sliding com- 
ponent represents only about one-haK of the gear losses, excluding bearing losses. The 
rolling and windage losses account for the other 50 percent and must be properly ac- 
counted for to obtain an accurate estimate of gear power loss. In addition, figure 8 
clearly shows that using the coefficient of friction to predict part-load gear efficiency is 
grossly inadequate, since the sliding power loss decreases to zero at low loads. 

Finally, figure 8 illustrates the fact that the rolling-element support -bearing losses 
are generally a significant portion of the total system losses and should not be ignored. 

In the case of the test gears that were analyzed, the bearing power loss varied from 
about 10 percent of the total system power loss at the lowest speed and highest torque to 
about 40 percent of the total at the highest speed and zero load. 

Figure 9 shows the power loss of the spur-gear system described in reference 6 as 
a percentage of the full-load power loss at several levels of loading. The lowest loading 
curve (an input torque of 0. 54 N-m (0. 4 ft-lbf)) represents the no-load or tare-loss 
curve. It is instructive to note that the tare power loss of an unloaded gearset theoreti- 
cally can reach 65 percent of the full-load loss (at an input torque of 271 N-m (200 ft-lbf) 
and a speed of 2000 rpm). This effect is due to the fact that at higher speeds the rolling, 
windage, and support -bearing losses are relatively insensitive to load but make up a 
major part of the system loss. As the load is decreased, only the load -dependent slid- 
ing loss decreases, and the system loss remains high. 


SUMMARY OF RESULTS 

A simple method for predicting the power loss and efficiency of a steel spur gearset 
of arbitrary geometry supported by ball bearings was developed. The method algebra- 
ically accounts for losses due to gear sliding, rolling, and windage and incorporates an 
expression for ball-bearing power loss. This method provides an accurate estimate of 
spur -gear -system efficiency at part load as well as full load. The analysis was com- 
pared with test data generated on a recirculating -power spur-gear test rig from another 
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investigation for a wide range of speeds and loads and for various methods of gear lubri- 
cation. The theoretical contributions of the individual component losses were compared 
with the total gear-system loss. The following results were obtained; 

1. The analysis generally showed good agreement with the test data at all condi- 
tions except low oil flow rates, where the analysis slightly underestimated the measured 
gear-system efficiency. 

2. The predicted sensitivity of power loss to change in speed and load closely re- 
sembled that determined experimentally. 

3. The contributions of gear rolling -traction loss and support-bearing loss to the 
total system loss were significant. The contribution of windage was also significant, 
but to a lesser degree. These losses should be accounted for in the loss calculation, 
particularly under part-load conditions. 

4. The unloaded or tare power loss of the gear system at operating speed could be 
as much as 65 percent of the maximum -load loss. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 10, 1979, 

505-04. 


17 



APPENDIX A 


EQUATIONS USED IN CALCULATIONS OF GEAR POWER LOSS 


Equations for Path of Contact 


The path of contact is a line in space mapped out by the point of contact of two gears 
in mesh. The path of contact is a straight line for involute gear geometry. The points 
Xj to X^ shown in figure 1 and Xp represent the following sequence of events; 

(1) Xj^ - start of mesh cycle, two teeth share the load 

(2) X 2 - start of single -tooth contact 

(3) Xg - end of single -tooth contact 

(4) X^ - end of mesh cycle 

(5) Xp - pitch point 

The following equations for these points are given in reference 3: 


g jp 

D 

P ^p 


\p=^p 


2C 


D = D + 

a, P P p 


10 


D = D cos 
o,g g 


2C 


D„ „ = D„ + 


10 


a, g g p 


h = 


D IT 

0>P 

N„ 


(D + D )sin e 

Y - F & 

^A 




y (Al) 


(continued) 
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/ 2 2 

^4=»-K<p-'’?,p) 

X2=^4-/'b 


1/2 


/ (Al) 


Xp=Xi + CjO 


'2 + N ^^2 ^ n ^^ 0\2 


ni/2 


i.f 


2./ 


2.f 


where 


Cj^Q = 0. 0254 (SI units) 

= 1.0 (U.S. customary units) 

Six additional lengths along the path of contact are required in the sliding and roll- 
ing power loss calculations: 


^1 = 


X1 + X2 


l2 = ^1 - Xj + X3 


f3=X4-Xs+X2-Xi 


f ^ - X2 + A-f 


^5=X3-^2 


f6=X4-Xi 


(A2) 


where 


A£ - 3xl0"^m(0.00l in.) 
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Sliding and Rolling Velocity 


The sliding velocity is the difference in surface velocity of the two gears at the 
point of contact. The sliding velocity varies along the path of contact and changes direc- 
tion at the pitch point. Since we are interested only in the magnitude of the sliding 
velocity, we eliminate the change of sign by taking the absolute value of X - Xp: 

Vg(X) = Vg - Vp = (cUg + cUp) |X - Xp I (A3) 

or in terms of pinion speed 

0. 1047(1 + mJn^lX - X^| 

Vg(X) = ^ ^ (A4) 

The rolling velocity as defined in the Benedict and Kelley friction coefficient (ref. 8) is 


V„ + V„ 


In reference 5, the rolling velocity is shown to be 

sin e - |X - Xp|(m - 1) 


VT=Vg.Vp=2V. 




This equation can be rewritten in terms of pinion speed as 


V.p(X) = 0.1047 npDp 


sin 0 - |X - Xp|(m - 1) 


D, 


(A5) 


Gear Load 

The normal load on the gear tooth surface may be calculated from the applied pinion 
torque: 
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2 


Vfr 


W. 


COS 6 


2T 


w 


n 


JL 


D cos B 
P ^ 


(A6) 


Between points X 2 and along the path of contact, w^^ is carried by one tooth. 
Between Xj^ and X 2 and between Xg and X^, the load is shared by two gears. Thus 


w(X) =w^ for Xg ^ X < Xg ^ 


w > 

w(X) =— 0 . forXj<X<Xg 

^ and Xg < X < X^ ^ 


(A7) 


Equivalent Rolling Radius 

The gear contact can be approximated by two equivalent rollers in contact. The 
equivalent roller radius changes along the path of contact and can be found as follows: 


D sin 6 ^ 


D sin 6 


R B^O. P 

eq R 




+ R 
o, g o,p 


(A8) 
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APPENDIX B 


FILM -THICKNESS EQUATION FOR CALCULATION OF ROLLING POWER LOSS 

Adaptation of Film Thickness Equation of Reference 11 to Gear Geometry 

The expression for central film thickness in an EHD contact developed in refer- 
ence 11 is 

= 2. 69 R^U®- 067(1 _ q. 61 e"°' (Bl) 

This general EHD equation is modified here to incorporate gear parameters by taking 
each nondimensional grouping and explaining its application to the gear contact: 

U speed parameter, uPq/E'R^ 

u surface velocity; in a gear contact, the average rolling velocity is used for u; 
since Vr^ is defined as V + V here, u = Vrp/2 

O P A 

n n 

/Iq atmospheric absolute viscosity of lubricant, N sec/m*^ (Ibf sec/in ) 

; if we assume that both the gear and the pinion are made of 

steel, E' = 2.276x10^^ N/m^ (33x10® Ibf/in^) 

equivalent rolling radius as defined in equation (A8), R (X) 

eq 

G material parameter, E'a, 4617 

a pressure-viscosity coefficient; for a typical mineral-oil gear lubricant, 
2.030x10"® m^/N (1. 4x10"'^ in^/lbf) 

W load parameter, Fj^/E'R^ 

Fj^ normal load on EHD contact; in gear contact, Fj^ =w(X), as defined in appendix A 
(under Gear Load) 

k elliplicity parameter, a/b; a and b define size of contact area; for gears in line 
contact, a constant value of 12 is assumed 

Combining these constants and variables gives 

h(X) =CjjlV,j,(X)Po|®-®’^|w(X)j-®-®®'^[R^q(X)J®-4®^ (B2) 

where 
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= 2.051x10"’^ (SI units) 

= 4.34x10'^ (U. S. customary units) 


Thermal Reduction Factor 


At high pitch line velocities, the film thickness calculated by equation (B2) over- 
estimates the actual film thickness. In reference 12, a thermal reduction factor is 
developed to correct for the inlet oil shear heating that occurs in high-speed EHD con- 
tacts. The thermal reduction factor (p^ is reproduced in figure 3 as a function of 
The value of can be calculated as follows: 


Q (X) =- 


' 12^0 


u“6 


K, 


(B3) 


where 

Cj2 = 1x10"^ (SI units) 

= 4. 629 (U. S. customary units) 

Since u is defined as the average surface velocity in reference 12, we must use 
u = Vrp(X)/2. The thermal conductivity of the lubricant is assumed to be constant. 
A typical value of Kj. for a mineral -oil lubricant is 

Kj = 0. 125 W/m K (0. 0725 Btu/hr ft °F) 

If the viscosity-temperature characteristics are known, 6 can be found for a given oil. 
The following equation must be solved for 6 by using known values of p at tempera- 
ture i3 and of Pq at temperature /3 q: 

At atmospheric pressure, p = 0; thus 


In Po - In p 
6 = ^ 

a -do 


For a typical mineral-oil lubricant. 
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Combining these values in equation (B3) yields 




where 


Cj3 =6.8x10“^ (SI units) 

= 0. 303x10"^ (U.S. customary units) 

After has been calculated, can be found in figure 3. 



APPENDIX C 


WINDAGE POWER LOSS EQUATION 


An expression for gear windage loss is developed here from data obtained in the 
studies of references 13 and 14 on turbine rotor windage losses. The moment coeffi- 
cient for a rotating disk is commonly expressed as 


m 


_C^ 

ipcV 

2 


(Cl) 


where 


C = 1.0 (SI units) 

0 

= 1.5x10“ (U. S. customary units) 

This coefficient was evaluated experimentally and found to be of the form 


C =— (C2) 

Re^ 

(refs. 13 and 14). This expression was modified to account for disk width as follows: 


C 


m 



(C3) 


where C^ q is the moment coefficient for a disk of zero thickness (ref. 13). This 
equation was verified for ratios of disk thickness to radius up to 0. 15. 

Reference 13 gives equation (C2) as 


'm, 0 


0.09 


Re 


0.21 


(C4) 


for the unbladed disk. Reference 14 gives equation (C2) as 

/S ^ 

0 . 102 ( 1 ') 

c. 


'm, 0 


Re 


0.20 


(C5) 
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where S is the axial clearance between the rotor and the housing wall, and R is the 
radius of the disk. (In ref. 14, the constant 0. 102 is printed as 0. 0102 in the expres- 
sion for C . Inspection of the data indicates that this constant must be 0. 102. ) Equa- 
tion (C5) applies when turbulent flow exists within the gearbox and the wall and gear 
boundary layers are not merged. 

In a gearbox application, generally there is significant axial clearance between the 
gear and the housing wall. For this analysis, S/R 0.217 (the largest value in ref. 14) 
is used. Thus equation (C5) becomes 


'm, 0 


0.0876 


Re 


0.20 


(C6) 


This value agrees with equation (C4) from reference 13. Combining equations (Cl), 
(C3), and (C6) results in 


m 


. ^ 14 ^ 

1 pc.2r5 

2 


= (1 + 2.3-t- 
R 


t\ 0.09 


Re 


0.2 


The drag torque T is 


T 


3J.yO.09 


R/r^O.2 


Converting this equation to a power loss and simplifying yield 


'^W =C4(l + 2.3-^y®- 


(C7) 


(C8) 


where 

C^ = 2.04x10'^ (SI units) 

2. 70x10"^^ (U. S. customary units) 
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TABLE I. - GEAR GEOMETRY AND OPERATING PARAMETERS 
FOR TEST GEARS OF REFERENCE 6 


Pitch diameter, cm (in.) 

Pinion 15.2(6) 

Gear 25.4(10) 

Number of teeth 

Pinion 48 

Gear 80 

Diametral pitch 8 

Pressure angle, deg 20 

Width, cm (in.) 4.0(1. 563) 

or 7.94(3. 125) 

Lubricant Mineral oil with 

antioxidant additive 

Viscosity at oil jet 

temperature of 333 K (140° F), cm^/sec 0.60 


TABLE II. - CONSTANTS USED IN GEAR POWER 


LOSS EQUATIONS 
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TABLE m. - SUMMARY OF POWER LOSS EQUATIONS FOR SPUR-GEAR SYSTEM AND DESIGN EXAMPLE 


[Gear data: N_, 48; N„, 80; S’, 8; 6, 20°; nip., 1.666; 0.0397 m (1.5625 in.); operating conditions; n„, 2000 rpm; T , 

P 6 c .2\ o V p 


271 N-m (2400 in. -Ibf); Mq, 0.05 N sec/m^ (7.25x10'” Ibf sec/in^); I'g, 0.60 cm'^/sec (6. 459x10”“^ ft^/sec); f^, 2; 
bearing data: D^, 0. 07 m (2. 75 in. ); C^, 17 436 N (3920 Ibf); Fg.p = w. 


(a) Equations for path of contact, windage loss, and bearing loss 


Parameter 

Symbol 

Formula 

Result 

Gear pitch diameter, m (in.) 

D 

g 

CioV^ 

0.2540(10.00) 

Pinion pitch diameter, m (in. ) 


CioV^ 

0. 1524(6.000) 

Pinion base circle diameter, m (in. ) 

®o,p 

D cos 6 

0. 1432(5. 638) 

Pinion tip diameter, m (in. ) 

D 

a,P 

Dp+2C,o/.y> 

0. 1588(6.250) 

Gear base circle diameter, m (in. ) 

®o, g 

D cos 0 
g 

0.2387(9.3969) 

Gear tip diameter, m (in.) 

^a, g 


0.2604(10.25) 

Base pitch, m (in. ) 

Pb 


9.372xl0'^(0.3690) 

Length between interference points, m (in. ) 

^A 

(Dp + Dg)sin 0/2 

69.50xl0'^(2.736) 

Start of double-tooth contact, m (in. ) 

Xl 

/ 2 2 

XA-O-H^L'^g) 

17.47xl0'^(0.6890) 

End of single-tooth contac, m (in. ) 

^3 

^1 V'b 

26. 84xl0‘®(1.058) 

End of mesh cycle, m (in. ) 

^4 


34.32xl0'^(1.349) 

End of double-tooth contact, m (in. ) 

^2 

^4-/'b 

24. 95xl0"^(0.9799) 

Pitch point, m (in.) 

Xp 

Xi + Cioj[(2 + Npmg)/a/]2 

26.03xl0'^(1.026) 



51 1/2 

- (Npmg cos 6nSr\ - CioNpmg sin 0/2./ 






Midpoint between Xj^ and Xg, m(in.) 


(X^ + X^)/2 

21.20xl0‘^(0. 8345) 

Midpoint between Xg and X^, m (in. ) 

h. 

fl-X^+Xg 

30.60xl0“^(1.204) 

Length of single-tooth contact, m (in.) 

^3 

X 4 -X 3 + X 2 -Xi 

14.96xl0'^(0.5818) 

Start of double -tooth contact, m (in.) 

^4 

X 2 + Af 

24.98xl0‘^(0.9809) 

Length of double-tooth contact, m (in. ) 

^5 

X3-X2 

18.94xl0‘^(0.0744) 

Total length of contact, m (in. ) 

^6 

X 4 -X 1 

16. 85xl0"^(0. 6599) 

Gear tooth normal load, N(lbf) 

^n 

2 Tp/(Dp cos e) 

3785(850.9) 

Gear windage loss, kW (hp) 

^W,g 

Cg(l + 2.3.>/Rg)(np/mg)2-2R4- 6 ( 0.028 

0.0164(0.0220) 

Pinion windage loss, kW (hp) 

^W,p 

Cg(l + 2.3 .¥/n^)n\- ^R^- ®(0. 028 Mg + Cg)®-^ 

0.0084(0.0112) 

Gear bearing load-dependent torque, N-m 
(in-lb) 

L, g 

0. 0009 

0.0351(0.3107) 

Pinion bearing load-dependent torque, N-m 
(in-lb) 

Mx 

L,p 

0.0009 

0.0351(0.3107) 

> 

Gear bearing viscous torque, N-m (in-lbf) 

M„ 

V,g 


0. 1079(1.024) 

Pinion bearing viscous torque, N-m (in-lbf) 

M.x 

V,p 


0. 1634(1. 437) 

Total gear bearing torque loss, N-m (in-lbf) 

M 

g 

^L,g-^V,g 

0. 1430(1.335) 

Total pinion bearing torque loss, N-m 
(in-lbf) 

M 

P 

^L,p"^V,p 

0.1985(1.748) 

Total support -bearing power loss, kW (hp) 

^BRG 

2 Cg(Mgng + Mpnp) 

0. 1194(0. 1620) 




TABLE m. - Concluded. 


CO 

to 


(b) Mesh losses along path of contact 


Parameter 

Symbol 

Formula 

Result at - 




X = fj 

— 

X 

II 

II 

X 

X 

II 

X 

Sliding velocity, 
m/sec (in/sec) 

Vs 

0.0147(1 + mg)np(X - Xp)/mg 

1.608(64. 17) 

1. 542(59 . 48) 

0.3418(15. 11) 

0 

Rolling velocity, 
m/sec (in/sec) 

Vt 

0. 1047 npDp[sin 6 - |X - Xp|(mg - 1)/Dg] 

10.51(413. 7) 

10. 53(414.9) 

10. 83(426.2) 

10.91(429.7) 

Tooth normal 
load, N (Ibf) 

Friction coeffi- 

w 

/ 

Wn or w/2 
(a) 

0.0127 log[CjW/(;j7-MQVgV2)] 

(1892)(425. 7) 
0.0279 

1892(425. 7) 
0.0282 

3785(851.3) 

0.0394 

3785(851.3) 


cient 







Sliding power 


CgVg/'w 

0.0852(0. 1155) 

0.0823(0. 1082) 

0.0420(0.0768) 


^s 


loss, kW (hp) 







Gear equivalent 

^.g 

Dp(sin e)/2+ |X-Xp| 

0.0309(1.218) 

0.0307(1.204) 

0.0271(1.071) 

0.0261(1.026) 

radius, m (in. ) 






Pinion equivalent 


Dg(sinfi)/2 - lx -Xpl 

0.0386(1.519) 

0.0388(1. 533) 

0.0424(1. 665) 

0.0434(1, 710) 

radius, m (in. ) 






Equivalent con- 
tact radius, m 


R R /(R + R ) 

o,g o,p^^ o,g o, p' 

0.0172(0.6758) 

0.0172(0.6742) 


0.0163(0.6413) 

% 


(in.) 

Film thickness, 
m (in. ) 

r fV u ■,0.67 -0.067p0. 464 
Cji(V.pM0> % 

1.247xl0'®(49. 20x10’®) 

1.248xl0’®(49. 24x10’®) 

1. 192xl0’®(47. 02x10’®) 

h 

Qm 


Thermal loading 

^la'^o^T 

0.376 

0.377 

0.405 


factor 

Thermal reduc- 

^t 

(b) 

0.931 

0.931 

0.924 


tion factor 
Rolling power 

^2^3VT*’‘^t'^ 

0.0436(0.0583) 

0.0437(0.0585) 

0.0429(0.0575) 


^R 


loss, kW (hp) 






®See fig. 1. 
fitr. 3. 



(c) Average mesh loss and system efficiency 


Parameter 

Symbol 

Formula 

Result 

Sliding power loss, kW (hp) 
Rolling power loss, kW (hp) 

Total system power loss, kW (hp) 
Power into gearset, kW (hp) 
System efficiency, percent 

Ps 

Pr 

p 

^TOT 

Pin 
. ( 

{|Ps(fi) + Ps(f2)J% + (Ps(^4»(%)/2} A 
Ar<V + Pr( VI^3 " (Pr(^P»%)A 
Ps Pr P\V, g P\V, p PrRG 

P P 

Pm ' Ptot/Pin 

0. 1511(0.2016) 
0.0823(0. 1095) 
0. 3776(0. 5063) 
56.79(76.2) 
99.34 








Power loss due to gears, hp 






50 100 150 200 

Pinion torque, ft-ibf 


Figure 5. - Comparison of predicted gearbox bearing 
power ioss with data of reference 6. 


CO 


Overali gearbox efficiency, percent 




(b) Pinion speed, 1000 rpm. 



50 100 150 200 


Pinion torque, ft-lbf 

(c) Pinion speed, 2000 rpm. 

Figure 6. - Comparison of predicted werall gearbox efficiency 
with range of data from reference 6. Gear width, 4.0 cen- 
timeters (1.563 in. ), 
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Figure 8. - Theoretical percentage comparison of sources of gearbox power loss as func- 
tion of input torque for three pinion speeds. Gear width, 40 centimeters (L563 in. I; 
gear geometry and operating parameters given in table I. 
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Figure 9. - Theoretical percentage of spur-gear-system power ioss at 
part-ioad conditions relative to full-load power loss. Gear width, 

4.0 centimeters (L563 in.); gear geometry and operating parameters 
given in table I. 
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